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Abstract 

The standard intensity-based approach for modeling defaults is generalized by making 
the deterministic term structure of the survival probability stochastic via a common 
jump process. The survival copula of the vector of default times is derived and it is 
shown to be explicit and of the functional form as dealt with in the work of Sibuya. 
t-H Besides the parameters of the jump process, the marginal survival functions of the 

default times appear in the copula. Sibuya copulas therefore allow for functional 
parameters and asymmetries. Due to the jump process in the construction, they 
allow for a singular component. Depending on the parameters, they may also be 
extreme-value copulas or Levy-frailty copulas. Further, Sibuya copulas are easy 
to sample in any dimension. Properties of Sibuya copulas including positive lower 
orthant dependence, tail dependence, and extremal dependence are investigated. An 
t-H application to pricing first-to-default contracts is outlined and further generalizations 

of this copula class are addressed. 
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1 Introduction 

(N 

{Sj A (i-dimensional copula is a (i-dimensional distribution function with standard uniform 

univariate margins. The main goal of the present work is to construct a flexible class of 
d-dimensional copulas based on a multivariate default model and investigate its properties. 
The multivariate default model considered is a generalization of the standard intensity- 
• • based approach by using a common jump process. After presenting this model, the 

survival copula combining the default times is explicitly derived. It is of the functional 
form as appearing in Sibuya (1959). The resulting class of "Sibuya copulas" is remarkable 



in many ways. A first distinguishing feature compared to other copula classes is the 
fact that Sibuya copulas have functional parameters and therefore quite flexible in terms 
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2 A default model and its implied dependence structure 



of its properties. As a second important feature, Sibuya copulas are not restricted 
to functional symmetry, i.e., exchangeability This is important in large dimensions 
when exchangeability becomes a more and more restrictive assumption. Third, due to 
the construction with a common jump process, Sibuya copulas may have a singular 
component. This is important for applications such as multivariate default models since 
it translates to a positive probability that several components of a system default at the 
same time. Another important feature of large-dimensional copulas in applications is 
sampling. Due to the construction of Sibuya copulas via a default model, it is easy to 
draw vectors of random variates from these copulas. 

Due to their general functional form, it seems difficult to determine the properties 
of Sibuya copulas in general. For investigating the properties, we therefore consider a 
working example throughout the paper. Even for this example, Sibuya copulas are seen 
to allow for many interesting properties. For example, the considered example is an 
extreme- value copula and also, as a special case, a Levy-frailty copula. Further, the lower 
and upper extremal-dependence coefficients may be derived explicitly. Moreover, for the 
bivariate case, the working example is of Marshall- Olkin type. 

As an application, we consider the valuation of a financial derivative contract, known 
as first-to-default swaps. We derive the pricing equation in the proposed dependence 
framework, and obtain an analytical formula for the fair spread of such type of contracts. 
Finally, possible extensions of the construction principle for Sibuya copulas are given. 

The article is organized as follows. Section 2.1 recalls the standard intensity-based 
approach for default modeling. A generalization utilizing a common jump process is given 



in Section [2. 2| In Section 2.3 we compute the joint survival function of the default times. 
The corresponding Sibuya type copulas are then derived in Section 2A In Section [3] the 
properties of this distributional class are investigated, including positive lower orthant 
dependence, tail dependence, extremal dependence, and sampling. Section [4] outlines an 
application to the pricing of first-to-default contracts, and Section [5] briefly addresses 
possible extensions of the construction principle for Sibuya copulas. Finally, Section [6] 
concludes. For the reader's convenience, proofs are given in the Appendix. 



2 A default model and its implied dependence structure 

In this section, we derive a dependence structure based on a default model in a stochastic 
intensity framework. First, the standard intensity-based model is reviewed. Next, the 
natural extension consisting of making the intensities stochastic is considered. The 
stochastic intensities are restricted to take the form of a non-decreasing deterministic part 
and a non-decreasing jump process. The latter, being common to all entities, generates 
the dependence among the constituents. 

2.1 The standard intensity-based default 

In the standard intensity-based approach, default times are modeled as the first-jump 
times of a (possibly non-homogeneous) Poisson process. In other words, the default time 
of the ith of d components in a portfolio is modeled by a deterministic, non-negative 
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intensity Aj : [0, oo) — > [0, oo), % G {l,...,d}. Given the intensity Aj, the survival 
probability pi(t) of component i until time t is given by 



pi(t) := exp(-Ai(t)), t G [0,oo), 



(1) 



for the integrated rate function Aj(t) := Jjj Aj(s) ds, t G [0, oo). The canonical construction 
of the default time fj is then given by 



f; :=inf{t > 0:pi(t) < U { }, 



(2) 



for triggers Uj ~ U [0, 1], i e see Bielecki and Rutkowski (2002 pp. 227) or 



Schonbucher (2003 p. 122). The survival function Si for the ith component at time t 



can now be computed as 

§i{t) := F(n >t)= ¥{pi{t) > Ui) = Pi{t), t G [0, oo). 

Usually, one seeks for generating default times that are mutually dependent. This 
is already introduced at various points in the literature, including, e.g., iLil (2000), 



Schonbucher and Schubert (2001 ), or Hofert and Scherer (2010) who introduce dependence 
by assuming a joint model for the vector of trigger variables (Ui, . . . , Ud) T - The following 
section we take another approach. We assume the triggers to be independent (see Section 
[5] for possible extensions) and introduce dependence by a common jump process. 



2.2 A generalized default model 

In what follows we generalize Mechanism Q for modeling the term structure of survival 
probability of an entity by considering Cox processes instead of Poisson processes. The 
deterministic survival process (pi(t))te[0,oo) ^ s replaced by a stochastic one, i.e., 

Pi {t) :=exp(-A M ), t G [0,oo), (3) 

where (-2Q,t)te[o,oo)> i G {1, . . . , d}, are right-continuous, increasing stochastic processes 
with independent increments and X^q = 0, % G {1, . . . , d} that we naturally refer to as 
the integrated intensity process ("IIP"). The idea underlying the stochastic extension Q 
of ([!]) is that samples of the default time may be generated by a more general model 
without altering the distribution. 

As a realistic model for the processes (-Xi,t)te [0,00)5 ^ G {1, . . . , d}, we consider 

X i>t := Mi{t) + J t , tG[0,oo), ie{l,...,d}, (4) 

where Mi(t) := J * m(s) ds, t G [0, 00), with a deterministic function m : [0, 00) — > [0, 00), 
analogous to Aj before, and (Jt)te[o,oo) is a right-continuous, increasing jump process with 
independent increments and Jq = 0. 

Remark 2.1 

The reader may find quite restrictive to focus on IIP of the form given by Q. However, 
it is a rather general form, as we now explain. In order for X^t to be an IIP, it needs to 
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be almost surely non-decreasing, otherwise the process Pi(t) may not be a proper survival 
process. Therefore, if we restrict Xn to be right-continuous, with independent increments 
and stationary, then by definition Xa is a Levy subordinator, which in turns, implies 
that any IIP satisfying these constraints are necessarily of the form jiit + Jt where fii 



is a non-negative constant, see, e.g., Cont and Tankov (2004 ( p. 88). The class of IIP 
we consider in Q is even more general than Levy subordinators in the sense that they 
can be non-stationary via the (possibly non- linear) functions Mj(i). This shows that the 
class of IIP defined by Q already covers an important part of the admissible processes. 

The intuition behind Q is that the individual term structure of the survival probability, 
modeled via Mi, is hit by the jump process {Jt)te[o,oo) which models common shocks 
affecting the components. The default of an entity i is then modeled similar to Q via 



n := M{t > : Pi {t) < Ui} 



(5) 



for Ui ~ U[0, 1], i £ {1, . . . , d}, where U is assumed to be independent of (Xi : t)te[o,oo)i 

ie{i,...,d}. 

Note that the default times Tj, i G {1, . . . ,d}, are naturally dependent due to the 
common jump process (Jt)t£\o,oo)- Our main goal is to investigate this dependence. We 
first focus on the case where U ~ U[0, l] d , i.e., the dependence is solely induced by 
(Jt)t£[o,oo)- An extension to nested or hierarchical dependence structures, as well as 
dependent triggers is discussed in Sections 5^ and |5.2| respectively. 

Since the dependence structure resulting from Construction ^ is quite general, we 
consider the following working example throughout this article. 

Example 2.2 (Working example) 

As a working example ("E"), consider (Jt)t&[o.oo) to be of the form 

Jt = HN t , t G [0, oo), 

for a non-homogeneous Poisson process (Nt)te\o,oo) with Nq = and N t ~ Poi(A(t)), 
where A(t) := Jq \(s) ds for a deterministic, non-negative function A : [0,oo) — > [0,oo) 
such that A(0) = lim^o = 0, and where H > is a constant. This choice corresponds 



to the model of Hull and White (2008) where the jump size is set constant and equal to 
H. 



2.3 The joint survival function 

In this section, we derive the joint survival function associated to our default model, 
namely S(t) := P(ti > ti,...,Td > td), t := (t±, . . . ,td) T £ [0,oo] d . An analytical 



expression of this survival distribution can be found in Vrins (2010) in the particular 



case defined by the working example. In what follows, we need the following lemma. 



Lemma 2.3 

Let ak £ [0, oo), k £ {0, . . . , d}, d G N, with a$ := 0. Further, let bk G 
such that bk+i = cbk for c£l and k G {1, . . . , d — 1}. 



k G {!,...,<*}, 
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(1) J2k=i a k = J2k=i(d — k + — a(fc-i)), where clua denotes the fcth order statistic 
of a k , k e {l,...,d}, i.e., a (1) < < a (d) . 

(2) Ea=i(! - & fe)( a (fc) - = (1 - c6d)a (d) + (c - 1) J2t=i hci(k)- 
Proof 



The proof is given in A.l □ 

The following theorem presents the joint survival function of the default times Tj, 
i £ {1, ...,d}. Here, ipy(x) denotes the Laplace-Stieltjes transform of the random 
variable Y at x, i.e., iby(x) := E[exp(— xY)], x £ [0,oo]. 

Theorem 2.4 

For i E {1, . . . , d}, let the survival process pi and the corresponding default time Tj be 
specified as in (§, @, and (|]), i.e., let 

n = inf{t > : ex V (-(Mi(t) + J t )) < U { }, i £ {1, . . . , d}. 

Then, the joint survival function of the default times is given by 

d ih J. _ Jf (d-i + 1) 

m-n-^p — *«.>. w 

where to := and Si(t) := P(Tj > i), t 6 [0,oo], i £ {l,...,d}, denote the marginal 
survival functions corresponding to S, given by Si(t) = exp(— Mi(t))ibj t (l), t £ [0,oo], 
i £ {l,...,d}. 

Proof 



The proof is given in A. 2 □ 

We may infer from (6) that the joint survival function S is of a form as dealt with 
in Anderson et al. (1992). It is parameterized by the corresponding marginal survival 



functions Sj, i £ {1, . . . ,d}, and the common jump process ( Jt)te[o,oo)- The following 
corollary presents the functional form of S for the case of our working example. 

Corollary 2.5 (Working example) 



Let us consider the setup of the working example, i.e., Example 2.2. In this case, the 
marginal survival functions are given by 

Si(t) ( = ] e W (-(Mi(t) +A(t)(l -e- H ))). 

Since the non-homogeneous Poisson process (Nt) t e[o,oo) nas increment distribution TVt — 
^*(t-i) ~ P°K-M%)) ~~ A(*(i-i))) an d since the Laplace-Stieltjes transform of Y ~ Poi(7) 
equals tpy(x) := exp(— 7(1 — exp(— x))), x £ [0, 00), 7 £ (0, 00), we obtain 

tf(t)®«p(-£(l - e-^ +1 ))(A(t (i) ) - A(^ 1} )) + (1 - e"*) £ A(f (i) )) 
^ i=l i=l ' 

d 

■nw (7) 

i=i 



5 



2 A default model and its implied dependence structure 



By applying Lemma 2.3||(1) with cij := A(tj), i G {1, . . . , d}, ([7]) can be simplified to 



= n - i + a (*«) - A (t(i-i)))Si(u), 



(8) 



where (p(x,y,z) := exp(—z(l — e~ xy — x(l — e~ y ))) , x,y,z£ [0, oo) denotes the jointure 
function as introduced in Vrins (2010); note that the function (1 — e~ xy — x(l — e~ y )) is 
non-positive for x G [1, oo), y G [0,oo) , and it is decreasing in y for any fixed x G [1, oo). 

with c := e H , en := A(U), and h := e - H ^- i+1 \ 



2.3 



Further, by applying Lemma 
i G {1, . . . , d}, the joint survival function in ([7]) can also be expressed as 

d x d 



5(t) ( ^exp((l- e ^)5:(l 



-H(d-i), 



= 1 



)A(t (i) ) nw 



(9) 



i=l 



2.4 The implied dependence structure 

With the joint survival function and the corresponding marginal survival functions at 
hand, one can derive the copula which provides a link between these two pieces of the 
multivariate default model. 

Corollary 2.6 

Let u := (m, . . . , Ud) T £ [0, l] d and let Si, i G {1, . . . , d}, be given as in Theorem 



2.4 



Further, let 5~ denote the generalized inverse corresponding to Si, i G {1, . . . , d}, and let 
Sr(u.)u\ denote the ith order statistic of S~(ui), i G {1, . . . , d}, i.e., Sr(u.)^ < • • • < 
Sr(u,)uy The copula C corresponding to the joint survival function (ml) is then given by 



C(u) = S(S;(u 1 ),...,S^u d )) = l[— 



, | ^J„-, > (1) 



(10) 



Since S i (t) < Sj (t), t G [0,oo), if and only if Si(t) < Sj(t), t G [0,oo), the diagonal 
corresponding to (10) is given by 



(■{« ") "''II : ! 

i=l 



S (i-D ( " 



Remark 2.7 

The copula C is in fact the survival copula of the vector of default times (n, . . . , t^) t . 
From the form above one recognizes that C is of the form of a class of distributions 
introduced by Sibuya (1959), who called 

d 



n 



s. (u.) (i) s. K)(j_i) 

^J s _ ( . (1) 

S. (u.) (i) 



the dependence function of the joint survival function S'(t). We therefore refer to class of 
copulas C as given in (10) as Sibuya copulas. 
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The idea of constructing a copula via a multivariate default model was recently applied 



by Mai and Scherer (2009a). In their work, (^Q,t)ie[o,oo) takes the form (A h .( t )) tg [ 0)OO ), 
i £ {1, . . . , d\, for a common Levy subordinator (A 4 ) 4g [ 0jOO ) combined with a rescaling 
of the time-clock via functions hi, i S {1, . . . , d}. This results in a non-stationary, non- 
decreasing stochastic process PQ,t)te[o,oo)- The rescalings hi are monotonically increasing 
entity-dependent functions derived from the riskiness of the entities, i.e., the subordinator 
time is passing more rapidly for riskier entities so that the survival process at some 
standard time point is lower for riskier entities than for safer ones. This approach also 
results in a tractable dependence model for defaults. However, the derived copula is 
restricted to functional symmetry, also known as exchangeability. This drawback is shared 
by many copula classes including the class of Archimedean copulas and also, partly, by 
nested Archimedean copulas. It becomes a more and more restrictive assumption in large 
dimensions since it implies that, e.g., all bivariate margins of the copula are equal. This, 
in turn, implies that, e.g., it is not possible to construct asymmetric (non-exchangeable) 
joint distributions if all margins are identical. Such restrictive properties are rarely 
observed, especially for large-dimensional portfolios. Sibuya copulas do not suffer from 
this drawback. 

Example 2.8 (Working example) 

From @ and @ we may infer that the Sibuya copula C in the setup of the working 
example is given by 

(F\ d d 

C(u) W n^-i+1, H, A(S7(u. ) w ) - A(S7 («.)«_!))) n «i 

i=l i=l 

= exp((l - e~ H ) ]T(1 - e- H ^)A(Sr(u) {t) )) f[ u t . 

^ i=l ' i=l 

If H = 0, then C(u) = H(u) := Ylf=i u i ^ s the independence copula. In the bivariate case, 
C can be written as 

C(ui,u 2 ) ( =exp((l - e~ H ) 2 A(S7{u.)^))uiu 2 

= exp((l - e- H ) 2 A(min{5r(m), S 2 (u 2 )}))u lU2 (11) 

From this functional form one may derive that C allows for a singular component, given 
by{^G[0,l] 2 |5r(n 1 ) = < S 2 -(n 2 )}. 

3 Properties of the copula 

In this section we investigate some properties of the copula C. We present results about 
positive lower orthant dependence, tail dependence, extremal dependence, and sampling. 



Due to the quite general form of a Sibuya copula C, see (10), it is difficult to investigate 



tail and extremal dependence, even under (E), i.e., for the case of the working example. 



In Section 3.2 and 3.3 we therefore work out the details under the additional assumption 
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("A"), which means Mj(i) = fat, i G {1, . . . , d}, and A(i) = At, t G [0, oo), where //j and 
A are non-negative constants. 

Let us first explore this case a bit. Under (A), the generalized inverses of the marginal 
survival functions Si, i G {1, . . . ,d}, are given explicitly by (u) = — log(ti)/Aj, where 
we define Aj := fa + A(l — e~ H ), i G {1, . . . , d}, for convenience. In this setup, a Sibuya 
copula can be expressed as 

(A) ^ A A A A 

CO) = J[uiip(d-i + l,H,log(u7 X/X ) {i) -log(nr A/A ) (i _ 1) ) 

n' 1 t -\/\.,(l-e- H )(l-e- H ( d -*1) , 1<y . 
ui(u. y (i) (12) 



/(i) 



with (u. A ^ A ')(o) := 1- The corresponding diagonal is given by 



d 



C{u, ...,u) ( =V l[<p(d-i+l, H, log u-^/^r 1 /^-!])) 

i=l 



= ^-Ad-e-^t^-e-^VX (13) 

i.e., a power function, where the subscript A[j] stands for the ith largest of Ai, . . . , A<2 
with A[q] := oo. 

The following remark addresses several properties of this copula. 

Remark 3.1(1) Assuming fa G (0, oo), i G {1, . . . , d}, if H = or A = 0, then C is the 
independence copula LT. Further, if H f oo and A f oo, then C becomes the upper 
Frechet bound copula M. We therefore conclude that Sibuya copulas allow to capture 
the full range of positive lower orthant dependence. 



(2) The Sibuya copula C as given in (12) allows for asymmetries and therefore more 
realistic dependence structures, especially in large dimensions. 

(3) Also note that this Sibuya copula is max-stable and therefore an extreme-value 



copula, see, e.g., Nelsen (2007 pp. 95), hence Sibuya copulas can be extreme- value 
copulas. 



(4) In the homogeneous case, i.e., fa =: fi for all i G {1, . . . , d}, ( 12 ) can be written as 



i=l M=l ' M=l 

_ A, l-c(l-e-"( d -<)) 
i=l 



where c := A(l — e H )/(fi + A(l — e H )), and the last equation follows by changing 
the order of multiplication in the second product of the second last equation. Thus, 



under (A), homogeneous Sibuya copulas are Levy- frailty copulas, see Mai and Scherer 
fl2009bp . 
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(5) In the bivariate case, (12) becomes 



■ ( 1—01 1— V21 

u\Ui = mrnluj U2,u\u 2 \ 



,1-02 ' 



where 



(1 



A/Aj S [0,1], i.e. a Marshall-Olkin copula, see 



Marshall 



and Olkin (1967). Note that in this case, one has explicit formulas for Spearman's 



rho and Kendall's tau, see, e.g., Embrechts et al. (2001), for the tail-dependence 



coefficients Nelsen (2007 p. 215), as well as for the probability of falling on the 



singular component, see, e.g., Nelsen (2007, p. 54). 



3.1 Positive lower orthant dependence 

It follows directly from Equation ( 



10) that C(u) > rij=i u ii i- e -> that C is positive lower 



orthant dependent, see Joe (1997 p. 21). For the bivariate case this property is also 
called positive quadrant dependence and it implies, by definition, that all measures of 
concordance such as Spearman's rho, Kendall's tau are greater than or equal to zero. By 



Equation (10), the dependence function directly controls this dependence since 



F{Ui<u 1 ,...,U d <u d ) 
P(C/i <«i)---P(C/<i<«d) 



n 

i=l 



(d-i+l) 



(i-i) 



(i) 



(1) 



> 1. 



(14) 



Since the left-hand side of ( 14 ) can be written as 



F(Uj < uj | U\ < ui, . . . , Uj-i < Uj-i, Uj+i < Uj+i, ... ,Ud<u d ) 

F(Uj <Uj) 

for any j £ {1, . . . , d}, one also says that the copula C has the "bad news propagation" 
effect. 



3.2 Tail dependence 

We pursue with a bivariate notion of association known as tail dependence. For Xi ~ Fi, 
i £ {1,2}, with joint copula C, the lower and upper tail- dependence coefficient \[ and 
X u , respectively, are given by 

A, := limP(X 2 < Fz(u) | X 1 < Ff («)) = lim C ^ u \ 

u.\,0 ulO U 



X u := limP(X 2 > F~{u) I Xi > F{{u)) = lim 



1 - 2u + C(u,u) 
1 - u 



where the limits are assumed to exist. By definition, the lower, respectively upper, tail- 
dependence coefficient tells us the likelihood, in the limit, that X\ and X 2 are both small, 
respectively large, simultaneously. Note that if (J7i, L^)" 1 " ~ C, then (1 — Ui, 1 — L^) 1 " ~ C, 
the survival copula corresponding to C. Therefore, the lower and upper tail-dependence 
coefficients interchange when going from C to its survival copula C. 
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In our case, C is the survival copula of the vector of default times (ti,T2) t . Thus, 
the lower, respectively upper, tail-dependence coefficient tells us the likelihood, in the 
limit, that the two default times t±,T2 are jointly large, respectively small. So upper 
tail dependence means that a joint default model with a Sibuya copula as dependence 
structure produces joint defaults within a short amount of time. 

If one assumes only the setup of the working example, then one can at least say that 

V=limuexp((l - e -lf ) 2 A(niiii{Sf (u), («)})). 

Thus, if at least one of the individual survival functions Si(t) is zero at some finite t, 
then A; = 0. Further, if A is bounded, then A^ = 0. Under (A), C is a Marshall- Olkin 
copula. Thus, the lower and upper tail-dependence coefficients are given by 

Xi = and A n = min{i?i, V 2 j = ? -, , W1 sr, 

max{/ii,/U2} + A(l - e- H ) 



respectively, see, e.g., Nelsen (2007 p. 215). It follows that H | or A | or maxjjiii, t 
00 implies X u 4- 0. So if we suppose, in the limit, that entity i E {1,2} is extremely safe, 
i.e., \i = fa + A(l - e~ H ) J, 0, then also X u J. 0. Further, if H f 00 and A f 00 or if H f 00 
and fj,i I 0, i € {1, 2}, then A u f 1. 

3.3 Extremal dependence 



The notion of extremal dependence was introduced by Frahm (2006). For Xi ~ F q 
i E {1, . . . ,ci}, with joint copula C, the lower and upper extremal- dependence coefficient 
Xi and A u , respectively, are given by 

e, : = limP(max{Fi(Aj)} < u\ min {^(Ai)} < u) 

uiO l<i<d l<i<d 



lim 



C(u, . . . , u) 



*4o 1 - C(l - u, . . . , 1 - u) 
P( mm {F,(A,)} > 1 

l<t<d 

C(l-u,...,l-u) 



e u : = limP( min {Fj(Aj)} > u I maxiFjfAi)} > u) 

w\l l<i<d l<i<d 



= lim 

ufi 1 — c7(w, ■ • • , it) 

where the limits are assumed to exist. By definition, the lower (upper) extremal- 
dependence coefficient tells us the likelihood, in the limit, that the largest (smallest) 
value of Fi(Xi), i E {1, . . . , d}, is small (large) given that the smallest (largest) value is. 
Applied to the setup where C is the survival copula of the default times, this means that 
the lower (upper) extremal-dependence coefficient tells us the likelihood that the smallest 
(largest) default time is large (small) given that the largest (smallest) is. Thus, given that 
the first default happened within a short amount of time, the upper extremal-dependence 
coefficient tells us the likelihood of all other defaults also happening within a short 
amount of time. The following proposition gives explicit formulas for ei and e u under 
(A). 
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Proposition 3.2 

Under (A), the lower and upper extremal-dependence coefficients are given by 



£l = 0, 



(A) 
£u = 



v , nm+1 \A ~ A(l - e- g ) S? = | f | +1 (l - e-^->)/A [t] 
ick.,d} d-A(l-e-^)Ef=i(l-e-^-))/A H ' 



respectively, where the sum extends over all 2 d subsets I of {1, ... , d}. 
Proof 



The proof is given in A. 3 □ 



3.4 Sampling 

The intuitive construction principle of Sibuya copulas via a default model can be used for 
simulation. Principally, we have to simulate the vector r := (n, . . . , r^) T of individual 
default times and then return (Si(ti), . . . , Sd(Td)) T , a vector of random variates from 
the Sibuya copula C. Sampling a vector r involves drawing a vector U ~ U[0, l] d and 
sampling a path of the jump process (Jt)te[o,T]i where T is such that p-iiT) < U{ for 
all i 6 {1, . . . , d}. Then, r is determined. Note that the number of occurrences to be 
sampled from the jump process depends on the given trigger variates t/j, i £ {1, . . . , d}, 
as well as on the deterministic functions Mj, i S {1, . . . , d}. The following algorithm 
describes the general sampling procedure of C. 

Algorithm 3.3 

(1) sample I7i ~ TJ[0,1], t € {1, d] 

(2) t hfi := 0, t := 0, k := 1, and I k := {1, . . . , 4 

(3) repeat { 

(4) sample the fcth occurence of the jump process (Jt)t£[o,oo) 

(5) find I C 4 : i e J ^ > Pi(t fe ) = exp(-(M i (t fe ) + J tk )) 

(6) for i g / { # find for all i e I 

(7) if (J7j < pi(ifc-) = exp(-(M i (t fc ) + Jt^J)) r» := t fc 

(8) else { # C/i G (pi{t k -) , Pi(t k ~i)) 

(9) find Tj on t fc ) via Tj := M^(- log C/j - J tfc _ 1 ) 

(10) } 

(11) } 

(12) := Ik\I # indices i for which Tj have not been determined yet 

(13) if (4+i = 0) break 

(14) else k := k + 1 

(15) } 

(16) return (5i(n), . . . , S d (r d )) T 
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3 Properties of the copula 



Under (E), (Jt)te[o,oo) i s given by Jt = HN t , t G [0, oo), for N t ~ Poi(A(i)). In this case 
we have Jt k = Hk, k £ {1,2, . . . }. Further, Step (4) of Algorithm 3.3 can be achieved 



with the following algorithm, see, e.g., Devroye (1986 p. 257). 
Algorithm 3.4 

(1) sample E ~ Exp(l) 

(2) tfik := tfik—i + E H 1 kth occurrence of a homogeneous Poisson 

(3) # process with unit intensity 

(4) tk '■= A~(th,k) # kth occurrence of a non-homogeneous Poisson 

(5) # process with integrated rate function A 



Note that under (A), Step (9) of Algorithm 3.3 boils down to setting r< := (— logUi — 
H(k — l))//ij and Step (4) of Algorithm 3.4 to := t/^fc/A. 

Example 3.5 

Let us consider the copula C as given in (11). Since under (A), the bivariate C is 
a Marshall- Olkin copula, we consider a more general example here. For this, let the 
"intensities" A and Hi, i £ {1, . . . , d}, be linear (instead of constant), i.e., let 

A(s) = a x s + b x , m(s) = aiS + b{, i £ {1, . . . , d}, 

where a\,b\,cii,bi S [0, oo), i S {1, . . . , d}. Further, let us assume the non-trivial case 
where not both oa (a«) and b\ (pi) are zero simultaneously. Letting c := (1 — e~ H ) we 
obtain 

A(t) = a x t 2 /2 + b x t, Mi(t) = ai t 2 /2 + ht, 

Pi (t) = exp(-( ai t 2 /2 + bit + HN t )), 

Si(t) = exp(-((aj + ca x )t 2 /2 + (h + cb x )t)) 

for i € { 1, . . . , d} . The corresponding inverses are 



A"(t) = t/b x l {ax=0} + (Vb 2 + 2at - b)/at {ax>0} , 
Mr(t) = t/bit {ai=0} + (Vb 2 + 2at - b)/al {ai>0} , 

Q - ( x _ V(bi + cb x ) 2 - 2(ai + ca x ) log(u) - (6j + c6 A ) 
\ u ) — 



Note that if tk denotes the feth occurrence of the non-homogeneous Poisson (A r f ) te [ 0)OO ) 
process with integrated rate function A, then 

Pi(t k ) = exp(-(a4/2 + bit k + Hk)), 
P7(tk~) = exp(-(a4/2 + 6^ fc + - 1))). 
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4 An application to the pricing of first-to-default contracts 



With these quantities one can evaluate the copula C and apply Algorithm |3,3| for drawing 
vectors of random variates from C. Figures [T] and [2] show two examples. Note that 
although we choose linear intensities, the resulting structures are quite different, e.g., see 
the shape of the singular components. Further, one may again infer that these copulas 
are able to capture highly asymmetric structures. 



C{u\, ui 





0.0 0.2 



0.0 0.0 



).4 0.6 0.8 

Ui 



Figure 1 Sibuya Copula (11) and singular component for Example (3.5) with H = 10, 
a x = 0.1, b x = 4, a = (ai,a 2 ) T = (1,100) T , and b = {b 1: b^y = (5,0) T (left). 
1 000 generated vectors of random variates from this copula (right). The lower 
and upper tail-dependence coefficients are given by A; = and \ u = 0.44, 
respectively. 



4 An application to the pricing of first-to-default contracts 

In order to illustrate how a joint model with Sibuya dependence structure can be applied, 
we consider a financial application, consisting in the valuation of contracts known as 
fir st-to- default (FTD) swaps. The key element of this contract is a basket, i.e. a pool 
of d > 1 entities. The first party, named protection seller agrees to pay the contractual 
counterparty known as protection buyer a fraction of the contract notional if the first 
default in the basket happens before the maturity of the contract. The fraction, named 
as loss given default, is defined by one minus the recovery rate of the defaulted entity. 
In order to enjoy that protection, the protection buyer is willing to pay the protection 
seller a given premium at prespecified payment dates (usually quarterly IMM), up to 
the maturity of the deal or the first default time, whichever comes first. After the first 
default in the basket or at the maturity, the latest, the contract stops. 

Mathematically, risk-neutral pricing theory proves that the above two payment legs 
are, per unit notional, equal to E[exp(— rr^){l — R^)t^ T <T y] (the default leg) and 
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C(Ui,« 2 





0.0 0.2 



0.0 0.0 



).4 0.6 0.8 

V, 



Figure 2 Sibuya Copula (11) and singular component for Example (3.5) with H = 10, 
a x = 0.1, 6 A = 4, a = (ai,a 2 ) T = (1,100) T , and b = {b 1: b^y = (0,0) T (left). 
1000 generated vectors of random variates from this copula (right). The 
lower and upper tail-dependence coefficients are given by A; = and X u = 1, 
respectively. 



sK[Jq exp(— rt)l{ T >t y dt] (the premium leg; for simplicity, we assume continuous pre- 
mium payments), where 7~(i) := minjri, . . . , r^} with distribution function F^., Ri 
denotes the recovery rate of the ith entity (with index "(1)" denoting the recovery rate of 
the first defaulting entity in the basket), s is the agreed spread (running annual premium), 
r is the (assumed constant) risk-free interest rate used for discounting the cashflows, and 
T is the contract's maturity in years. A basic calculation leads to the present value of 
such a contract from the protection buyer's perspective, given by 



FTD T = E 



exp(-rr (1) )(l - < T} 



- sE 









exp(-rt)l { >t} dt 



= (1 - R (1) ) ^ T exp(-rt) dF T{1) (t) - s ^ exp(-rt)(l - F T(1) (t)) dt. 

Another product, known as credit default swaps ("CDS") allows one to derive risk- 
neutral default probability curves from market prices. This provides us with d univariate 
default distributions 1 — §i(t), t £ [0, oo], i e {1, ... ,d}. However, this is not enough 
to infer F T ^y the distribution of the first default time T(i). The survival distribution 
of rm is given by the diagonal of the joint survival distribution, which we cannot infer 
solely based on the d univariate statistics at our disposal. In order to fill that gap, one 
needs to assume a dependence structure between the default times, while preserving the 
(market implied) individual survival probabilities. Here we achieve this with the help of 
our derived Sibuya copula class. In the following, we shall see how one can derive the 
FTD default distribution. Interestingly, in our framework, the later will be shown to 
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4 An application to the pricing of first-to-default contracts 



admit a closed form expression; this is remarkable noting that this is not the case even 
for the "simple" Gaussian copula model. Since the purpose of this paper is not to discuss 
how the copula parameters can be calibrated and in order to avoid entering the technical 
details related to the application, we assume them as given (e.g., by an expert). 

For simplicity, assume the individual default times have the same survival function 
S(t) = exp(— At) calibrated to CDS market quotes. Technically, this corresponds either 
to a one-point or to a flat CDS spread curve. For the joint dependence structure of the 



vector r of default times, we use the Sibuya copula C as given in (12). Its diagonal is 



given in (13). Applying Integration by Parts, FTDy can be computed as 

FTD T = (1 - R (1) )(l - eM~rT)C(S(T), S(T))) 

- ((1 - R {1) )r + s) [ C(S(t), S(t)) exp(-ri) dt, 
Jo 

so that the fair spread, obtained by solving FTD^ = with respect to s =: s*, is given 

by 

_ t l-e-*-C(S(n...,S(T)) _ \ 
{) \f^C(S(t),...,S(t))exp(-rt)dt J 



Writing j3 := d — A(l — e H ) 2~Zf=i(l ~ e H<yd f° r the exponent appearing in (13), 

we have C(S(t), . . . , S(t)) = exp(—/3Xt), so that 

s* = (1 - R (1) )P~\. (15) 

Now let us consider the attainable FTD spreads. The boundary of the (H, A)-space, 



i.e., H = or A = 0, corresponds to independence, see (12). In this case (3 takes its 
largest value ((3 = d), so does the spread in ( J]~5| ) . The lowest attainable FTD spread 
is obtained for the smallest value of /3 {f3 = 1), which, in turn, is obtained by letting 
H f oo and A f oo. The corresponding limiting copula is seen to be the upper Frechet 
bound C{u) = M{u). This corresponds to maximum correlation where all entities default 
simultaneously when the first jump of the process (Jt)t&[o,oo) takes place. 
An example of FTD spread level curves is shown in Figure [3j 

Remark 4.1 

The joint default model which led to the Sibuya copula is built on a general default model, 
assuming stochastic default intensities of some entities. Those should not be confused with 
the deterministic intensities of the FTD basket constituents, bootstrapped from the CDS 
market. In other words, there is no obvious link between the functions S~ , % £ {1, . . . , d}, 
involved in the copula definition (and solely used to build the dependence function) and 
the univariate survival functions of the constituents of the basket, which will be plugged 
in the copula in order for the joint survival function to have the market-implied margins. 
However, in the particular case where one sets the model parameters S~ , i G {1, . . . , d}, 
such that the functions Si are equal to the (CDS-based) survival probability curves Si, 
i £ {1, . . . ,d}, then the link becomes obvious, namely each realization of the stochastic 
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4 An application to the pricing of first-to-default contracts 




Figure 3 We consider a basket consisting of five entities with individual 
default intensities Aj, i £ {1,...,5}, given by 0.05, 0.133, 

0. 133, 0.142, and 0.145, respectively. Further, we take A 
to be the mean and ^ to be the minimum of these values, 

1. e., A = 0.1206, m = 0.05, i G {l,...,d}. We also take 
R = 40%. The dotted lines indicate the boundary of the 
(H, A)-space. Each point on the remaining curves represents 
the corresponding FTD spread level s* for the associated 
parameter values H and X(H, s*). The level curves are simple, 
and an optimization on such a space is nicely behaved. 
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5 Generalizations 



survival process Pi(t) can be seen as a possible default probability curve for the ith entity, 
to which some probability is assigned, and such that E[pj(t)] = Si(t). Thus, the intuitive 
framework based on random intensities translates to the basket constituents only in this 
case. 



5 Generalizations 

A key feature of our model is that it can be easily extended without affecting its tractability. 
In this section we briefly investigate some generalizations of the copula construction 
via the default model specified in ([3]), Q, and ^ and show that it remains perfectly 
workable. Although a lot of such extensions are possible (for instance, the scalar H 
could be made name-specific so that the jumps are rescaled in a name-specific fashion, 
or replacing the jump process by a sum of jump processes each with possibly different 
scaling coefficients, etc.), we restrict ourselves to analyze two possible generalizations. 
First, we replace the single jump process (Jt)t£[o,oo) by a hierarchy of jump processes 
( < ^i^)ie[o,oo)) j € {1, . . . , J}. Then, we address the case of dependent trigger variables U. 



5.1 Generalization to hierarchical jump processes 

The jump process (</t)te[o,oo) is fatal in the sense that it hits all components of the 
default model simultaneously. In practical applications, it might be the case that only 
certain subgroups of components get hit. Such a hierarchical or sectorial behavior may be 
modeled via (possibly dependent) jump processes (Jj,t)t£[o,oo) f° r J € {1, ... , J}, where 
J is the number of sectors or subgroups. These subgroups often arise naturally from the 
application considered, e.g., by given industry sectors, macroeconomic effects, geographical 
regions, political decisions, or consumer trends. A default model incorporating such 
hierarchies can be constructed with the stochastic processes (A"j ]i ) tg [ oo ), i G {1, . . . , d}, 
in Q being replaced by 

X, i)t := Mi(t) + J jtt , t€[0,oo), je{l,...,J}, ie{l,...,dj}. 

The corresponding individual survival processes are then given by Pji{t) := exp(— Xja), 
t G [0, oo], and the default time of entity i in sector j by Tj\ := inf{t > : pji(t) < Uji}, 
where Uji ~ U[0, 1], independent for all j G {1, . . . , J}, i G {1, . . . ,dj}. For simplicity, 
we assume the jump processes (Jj,t)te\o,oa)j 3 e {!> ■■ ■ ■> to be independent in what 
follows. Then the joint survival function S can be derived similarly as in the proof 



of Theorem 2.4 First note that the individual survival functions Sji are given by 
S 3i {t) = exp(-M i (i))V'J J , t (l), t e [0,oo], j G {1,...,J}, i G {1,...,^}. The joint 
survival function S can then be calculated as 

J •' . ./., _ (dj -i+ 1) 
S(tn, . . . ,h dl ,. . . ,tji, . . . ,tjdj) = Y]_ II ~ — > 3 1 q\ Sji(tji), 

j=li=l '^ J:i ' t j{iV ' 
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5 Generalizations 



where tju\ denotes the ith smallest value of all components {tji, . . . ,tjd-} in sector j. 
The corresponding copula C is thus given by 



^) = nn 

j=i i=i 



-J (dj - i + l, 

j,S~ (u . ) j,S. (u.) v J ' 



4>J. g - ( . (1) 

j,S . (u . ) 



Uji 



Ji- 



which is a product of Sibuya copulas and therefore itself of the Sibuya type. Therefore, 
Sibuya copulas are able to capture such hierarchical default dependencies. 



5.2 Generalization to dependent trigger variables 

Now let us introduce dependence among the default triggers via U ~ Cu, i.e., the trigger 
variables Ui, i G {1, . . . , d}, are dependent according to the copula Cjj. Note that this 



does not influence the marginal distributions Si, i G {1, . . . ,d}, as given in (16). Redoing 
the calculations as carried out in the proof of Theorem |2.4| leads to joint survival function 



S(t) = E[C [/ (exp(-(Mi(ti) + J h ), . . .,exp(-(M d (t d ) + J td )))] 
= E 



r , I Si(ti) S d (t d ) , 



^(1) 



The corresponding copula C is therefore given by 
C(u) 



E 



r , e M-J s - iul )) 

c u[ — ; \^ ui,' 



exp(— J 



S-l ("l) 



-u d 



Example 5.1 

Let us assume the homogeneous case, i.e., assume that, pointwise, = Mi for all 
i G {2, . . . , d}. This implies that, pointwise, Si = Si, i G {2, . . . , d}. As an example 
where the copula C is given explicitly, consider Cjj(u) = aM(u) + (1 — a)U(u) for 
a G [0, 1], i.e., Cu{u) is a convex combination of the upper Frechet bound copula M and 
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the independence copula II. In this case, applying Theorem 2.4 leads to 



S(t) = E[amin{exp(-(Mi(ti) + J u ))} + (1 - a) ]J exp(-(Mi(ti) + J u )) 



i=l 



a 

aE[min{exp(-(M 1 (t i ) + J t .))}] + (1 - a)E[j[exp(-(Mi(fi) + <4)) 



aE[exp(-(M 1 (^ ) ) + J t(d) ))] 

d V"J t -J, (d-i + 1) 

+ (i - a) n w — 



a 



<* Vj t - J t (d-i + 1) 

The copula corresponding to £ is therefore given by 

C(u) = aM{u) + (l-a)n ^ 



8=1 



^v, , (i) 



-j (d-i + i; 

aM (u) + (i - «) n Si ( " (d " + ^ )) Si ( " (d " + ( 2 }' 



where u^+i := 1. Thus, we recognize that C is a convex combination of the upper Frechet 
bound copula M and the Sibuya copula as given in ( 10 ) for the homogeneous case. 



6 Conclusion 

We introduced an intuitive default model which extends the classical intensity-based 
approach by allowing the survival process to jump downwards, i.e., to be stochastic. We 
then derived the survival distribution and, as a corollary, the associated copula which is 
proven to be of Sibuya type. For that reason, they are named Sibuya copulas. Since the 
parameters of the marginal survival functions of the default times appear in the copula, 
Sibuya copulas allow for asymmetries. Due to the jump process in the construction, they 
allow for a singular component. We also showed that Sibuya copulas may be extreme- value 
copulas or Levy-frailty copulas, depending on the functional parameters chosen. From 
the construction principle presented, a sampling algorithm for these copulas is derived. 
Further, properties including positive lower orthant dependence, tail dependence, and 
extremal dependence are investigated. A financial application consisting in the pricing of 
first-to-default swaps is given, and the nice-looking related expressions further emphasize 
again the interesting tractability of the model. This tractability most likely results from 
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A Proofs 



the relatively simple form of the integrated intensity process (although corresponding 
to a quite general setup) combined to the Sibuya property. Finally, we showed that 
the dependence model easily extends in various ways, and, as an illustration, explicitly 
derived the construction of Sibuya copulas for two of them. 

A Proofs 

A.l Proof of Lemma IZ31 



Proof 

For Part 
k)an.\ implies 



(1)1 £g=i(<* " k + 1)^) = Etl(d - k)a (k) = £fcftd - k)a (k) = T,k=l(d 



d d d 

J2(d -k+ l)(a (fc) - a (fe _i)) = ^(d -k + l)a (fc) - ^(d - k + l)a(&-i) 

k=l k=l k=l 

d d d 

= J2(d -k + l)a (fc) -J2(d- k)a {k) = a (k) 

k=l k=l k=l 

d 

fe=l 



For Part 
c T,k=i b k a (k) imply 



( 2 )[ Efc=i( a (fc)- a (fc-i)) = a (d) an d Efe=i ha(k-i) = Efc=o & fc+i a (fc) = Y?k=i b k+\a(k) 



/",(! _ h)(a(k) - a (k-i)) — a (d) - ^2 M a (fc) - a (k-i)) 

k=l k=l 



d d 
fc=l fc=l 

- bdU(d) ~ bka (k) +C Y1 b k a (k) 
k=l k=l 

d-l 

(1 - b d )a {d) + (c-l)J2 ha (k ) 

k=l 
d 

(l-cb d )a {d) + (c-l)^2b k a {k) . □ 

k=l 



A.2 Proof of Theorem 1231 



Proof 

First consider the ith marginal survival function S{(t) = P(rj > t), t £ [0, oo]. By 
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conditioning on Jt, it can be computed via 

Si(t) = F( Pi (t) > Ui) = P(exp(-(Mi(i) + J t )) > Ui) 

= EpP(exp(-(Mi(t) + J t )) > Ui | J t )\ = E[exp(-(M i (t) + J t ))\ 

= exp(-Mi(t))E[exp(-J t )] = exp(-Mi(t))^ Jt (l). (16) 

Given this result and Lemma 2.3 (1) with a, := J^, z G {1, . . . , d}, the joint survival 
function of the default times can be computed via 

S(t)=F(p i (t i )>U i , »€{l,...,d}) 

= P(exp(-(M i (t i ) + J ti )) >[/i,iG{l,..., d}) 

= EpP(exp(-(Mi(ti) + J*,)) >Ui, ie{l,...,d}\J ti , i€{l,..., d})] 

d x n d 



E 



E 



i=l 



exp 



H exp(-(Mi(ti) + J ti )) = E exp - £ J t . JJ exp(-M;(*i)) 



^(d-i+l)(J t - J t 



(i-i), 



8=1 



i=l 
d 



i=l 



n 



d ibj. _ Jt (d-i + 1) d ib Jt —j f (d-i + 1) 

n 1,1 — -*w = n "' r 11 ',, — W.>. 



where, in the second last equality, the independence assumption between any non- 
overlapping increments of the jump process is used. □ 



A. 3 Proof of Proposition |3.2 



Proof 

The survival copula C corresponding to C can be recovered from C by the Poincare- 
Sylvester sieve formula. This implies 



C(l-u, 



!"«)= E ( 

IC{l,...,d} 



-l) w C(u 



, u 



j2 (-i) izi n^ i/W (^ Ai/(o/A ')Uf - H ^- H{d -' ) \ 



IC{l,...,d} 



1=1 



where denotes the indicator of i being in /, i G {1, . . . , d}. Now consider the term 
( u -Al/(.)/A.) j£ tliere i s no A; G {1,. . .,d} such that k <£ I, i.e., if I = {1,. . . ,d}, then 

(«~ A1/ (")/ A -)(^ = (u~ A/A ) (i ) = u~ A/A W. If there is precisely one k G {1, . . . ,d} such that 
^ /, then (u _A1/( ' ) / A )(j) = u~ A / A M if i G {2, . . . , d} and 1 if % = 1. If there are precisely 
two different h,k 2 in {l,...,d} such that h,k 2 $ J then (n- Alj (-)/ A ) (i) = u~ A / A M if 
i G {3, . . . , d} and 1 if i G {1, 2}. Continuing this way one obtains that 

(n -AM-)A. } J*. 

V ' W \ U - A / A W, iG{|/| + l,...,d}. 
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This implies that 

60. i- u )= yi ri u-^w)^)^ 4 - ) 

7C{1,...,4 i=|I|+l 
/C{l,...,d} 

= £ (_ 1 )|/| u W- c Em +1 «i j 

IC{l,...,d} 

where c := A(l — e~ H ) and Oj := (1 — e~ H ( d ~ l >) / \u\, i G {1, . . . ,d}. Finally, note that 

C(u,...,u) =u d ~ c ^ti a \ 
With these two ingredients, we obtain 



v u 
ei = iim 



d—c^_ 1 <H 



IC{l,...,d} 



d 

l =|I| + i ai 



/C{l,...,d} 

e u = lim 



E (_i)|/|„l'l- c £t m+ i* 



1 1 — U *-ii=l 1 



The first statement now directly follows from the formula for ei. For the second statement, 
apply l'Hopital's Rule. □ 
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